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ABSTRACT 


A  combustion  model  and  burning  rate  equation  for  an  aluminized 
polybutadiene  ammonium  perchlorate  composite  solid  propellant  was  developed 
in  terms  of  the  rate  of  heat  feedback  to  the  propellant  surface.  It  was 
found  that  the  burning  rate-pressure  relationship  for  this  propellant  was 
well  represented  by  an  equation  of  the  form: 
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where  Tp  c  and  d  are  constants.  The  feedback  temperature  in  the  equation, 
T 2,  was  determined  as  a  function  of  pressure  using  propellant  ingredient 
physical  chemical  properties  and  subscale  burning  rate  data.  The  equation 
was  used  to  calculate  the  propellant  extinction  behavior  in  test  chambers 
with  low  (or  subsonic)  exhaust  flow,  and  the  C*  extinction  in  test  chambers 
exhausting  with  supersonic  flow  (i.e.  rocket  motors).  These  calculations 
agreed  well  with  the  observed  extinction  characteristics  of  tiro  propellant. 
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The  purpose  of  this  paper  is  to  present  a  model  and  burning  rate 
equation  developed  to  describe  the  pressure  effects  on  propellant  burning 
rate  and  the  propellant  extinction  characteristics  of  an  aluminized  poly¬ 
butadiene  ammonium  perchlorate  composite  propellant.  The  analysis  was 
performed  using  burning  rate  data  measured  by  the  Aerojet  General  Corpora¬ 
tion,  Sacramento,  California,  in  subscale  tests  on  cured  strands  about  one 
inch  in  diameter  (2  to  200  psia)  and  full  scale  tests  in  a  5000  pound  motor 
(6  to  500  psia). 

The  propellant  has  a  high  solids  loading  and  the  primary  fuel  is 
aluminum.  The  burning  rate  is  controlled  by  varying  the  quantities  of 
coarse  and  finely  ground  oxidizer  used  in  the  particular  batch  of  propellant. 
Using  many  large  scale  batches  of  the  same  propellant  formulation  burning 
rates  were  measured  in  full  scale  and  subscale  tests  and  the  results  were 
assessed  to  determine  reproducibility  and  control.  It  was  observed  that  the 
burning  rate  was  strongly  dependent  on  the  design  of  the  apparatus  in  which 
the  burning  rate  was  measured  (i.e.  strand  bomb,  50  lb  motor,  large  motor) 
under  conditions  where  propellant  erosion  was  not  an  important  factor.  In 
particular,  the  design  with  the  best  conditions  for  heat  feedback  to  the 
burning  surface  resulted  in  the  highest  burning  rate  at  a  given  pressure. 

A  burning  rate  equation  directly  treating  the  rate  of  heat  feedback  was 
developed  and  the  parameters  determined  from  the  accurately-known  burning 
rate  data. 
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This  work  differs  from  previous  work  on  propellant  burning  rates  in 
that  the  burning  rate  is  related  directly  to  the  radiation  and  conduction 
contributions  of  the  energy  fed  back  to  the  propellant  surface  from  the 
flame.  The  treatment  reported  in  the  literature  providing  the  best  des¬ 
cription  of  propellant  burning  rates  over  wide  pressure  ranges  is  due  to 
M.  Summerfield. 

The  granular  diffusion  model  of  solid  propellant  burning  developed  by 
Summerfield  and  his  associates1  which  relates  chemical  reaction  and  gas 
diffusion  effects  to  propellant  burning  rates  has  been  applied  to  correlate 
and  predict  the  burning  rates  of  many  ammonium  perchlorate  propellants , 
and  to  interpret  the  mechanism  of  ammonium  perchlorate  composite  propel- 
lants.  The  equation  describing  propellant  burning  at  high  pressure  based 
on  the  granular  diffusion  model  is: 

1  _  a  ,  b 

r  “  p  (2) 

where  a  and  b  are  constants  determined  by  fitting  the  observed  burning  rate 
data. 

In  applying  Summerfield' s  model  to  las'  pressure  (below  10  psia) 
propellant  burning,  the  effects  of  premixed  chemical  reaction  time  of 
ammonium  perchlorate  decomposition  products,  ammonia  and  perchloric  acid, 
are  taken  into  account. 

The  analysis  presented  here  is  based  on  the  hypothesis  that  the  burn¬ 
ing  rate  of  a  solid  propellant  is  determined  largely  by  the  thermal  proper¬ 
ties  oi  the  propellant  components,  and  that  the  chemical  reaction  and  gas 
diffusion  effects  are  not  rate  determining.  That  is,  the  chemical  reactions 
are  assumed  to  be  rapid  at  temperatures  where  tire  propellant  ingredients 
gasify,  or  as  the  gaseous  products  are  heated  by  the  extremely  high 
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temperature  gradients  in  the  flame.  The  model  described  below  was  developed 
using  as  inputs  the  results  of  thermochemical  calculations  of  propellant 
properties  (e.g.  flame  temperature,  C*,  etc.)  and  data  on  the  temperatures 
at  which  the  major  ingredients  in  the  propellant  under  consideration  gasify 
at  rates  comparable  with  the  burning  rate  as  follows: 

Gasification  or 

Major  Propellant  Ingredient:  Decomposition  Temperature 

o  Polybutadiene  binder..  , 

(and  plasticizer) '>>  500°C 

o  Ammonium  Perchlorate5  600°C 

o  Aluminum5  1000  -  1S00°C 

depending  on  the 
partial  pressure 

As  shown  below,  the  burning  rate  equation  developed  from  this 
model  describes  the  observed  burning  rates  at  pressures  from  below  1  psia 
to  1000  psia,  and  in  addition,  provides  a  method  for  relating  subscale  and 
full  scale  data,  and  predicts  propellant  extinction  properties. 

A  schematic  representation  of  the  combustion  model  used  in  tins 
work  is  shown  in  Figure  1.  The  general  form  of  this  model  for  solid  pro- 
pellant  burning  has  been  verified  by  numerous  experiments,  as  summarized  in 
Reference 

It  was  assumed  that  in  the  burning  process,  solid  particles  of 
ammonium  perchlorate  and  aluminum  are  released  from  the  surface,  and  the 
solid  surface  recedes  at  the  rate  of  binder  pyrolyzes.5  Due  to  the  thermal 
conductivity  and  heat  capacitance  properties  of  the  binder,  aluminum  and 
ammonium  perchlorate ,  the  heat  arriving  at  the  solid  surface  generates  gases 
initially  at  S00°C  from  the  binder  while  the  aluminum  and  ammonium  perchlor¬ 
ate  are  at  temperatures  closer  to  20i  *C  when  they  leave  the  surface. 


In  Zone  1,  the  pyrolysis  gases  from  the  binder  are  heated  to  about 
1200°C  (determined  by  fitting  the  data)  while  the  solid  aluminum  and  solid 
ammonium  perchlorate  are  heated  to  an  average  temperature  of  500°  -  600°C. 
At  600°C,  under  the  transient  heating  conditions  present  in  burning  pro¬ 
pellant,  the  ammonium  perchlorate  gasifies  rapidly,  thus  identifying  the 
end  of  Zone  1.  It  was  assumed  that  the  binder  pyrolysis  reaction  and  some 
endothermic  decomposition  of  the  binder  gases  are  the  only  reactions 
occurring  in  Zone  1  and  that  the  heat  entering  Zone  1  by  conduction  and 
radiation  from  the  combustion  products  is  all  transferred  to  the  propellant 
surface,  or  goes  to  decompose  the  binder  and  raise  the  temperature  of  the 
binder  products,  aluminum  and  ammonium  perchlorate. 

Since  Zone  1  is  defined  in  terms  of  the  temperature  of  tire 
ammonium  perchlorate,  the  thickness  of  Zone  1  depends  on  the  incident  heat 
flux,  the  pressure  and  the  mass  of  materials.  The  resulting  thickness  is 
essentially  the  path  length  necessary  to  conduct  the  heat  from  Zone  2  back 
to  the  surface.  There  is  little,  if  any,  contribution  to  the  heat  transfer 
in  Zone  1  from  rad.ation,  since  Zone  1  consists  largely  of  particles  of 
oxidizer  and  fuel  which  intercept  most  of  the  radiant  energy. 

Zone  2  starts  when  the  ammonium  perchlorate  begins  to  gasify 
rapidly  and  ends  at  the  temperature  from  which  heat  is  fed  back  to  the  solid 
surface.  At  high  pressure,  heat  is  fed  back  from  products  at  the  flame 
*  temperature  (which  is  lower  than  the  aluminum  boiling  point) .  At  inter¬ 
mediate  and  low  pressures,  heat  is  fed  back  from  the  region  in  which 
aluminum  boils  rapidly  (1)  in  equilibrium  with  aluminum  vapor  at  the  burn¬ 
ing  pressui'e.  or  (2)  at  a  partial  pressure  lower  than  the  burning  pressure. 
It  is  assumed  that  at  intermediate  and  low  pressures  the  aluminum  liquid  and 
vapor  form  a  layer  between  the  flame  zone  and  the  propellant  surface, 
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shown  in  Figure  1,  and  the  temperature  of  this  layer,  controls  the  rate  of 

heat  feedback  to  the  propellant  surface. 

In  Zone  2,  tire  oxidizer  gasifies  and  the  oxidizer  gases  and  binder 

pyrolysis  products  react.  It  is  assumed  that  the  thickness  of  Zone  2  is 

determined  primarily  from  the  volume  occupied  by  pockets  of  the  oxidizer 

gases  at  the  average  temperature  in  the  zone.  Tire  oxidizer  particles 

gasify  rapidly  once  they  reach  about  600°C  and  form  spheres  of  gas.  These 

3 

spheres  occupy  a  volume  of  where  is  the  diameter  and  is  related 
to  the  pressure  by  the  ideal  gas  law  (Equation  3) 


where  N  is  the  number  of  moles  of  gas  produced  by  the  oxidizer  particle, 
R  is  the  gas  constant,  T  the  temperature,  and  p  the  pressure. 

The  chemical  reactions  between  tire  pyrolyzed  binder  and  the 
gasified  oxidizer  are  assumed  to  be  rapid  at  the  temperatures  in  Zone  2, 
since  many  of  the  species  are  free  radicals  (e.g.,  H,  0  and  hydrocarbon 
radicals).  The  thickness  of  Zone  2  (J^)  is: 


where  1^  is  a  constant  (determined  in  practice  from  high  pressure  fullscalc 

motor  burning  rate  data)  and  T  is  the  average  temperature  in  Zone  2 

T  +  T 
Z  _  2  1 

T  = - 2 - . 

Tire  surface  recession  rate  (r  )  of  the  binder  in  the  propellant 

rU/ 

under  consideration  (equal  to  the  burning  rate,  r)  is  expressed  as  a 
function  of  the  heat  flux.,  incident  on  unit  area  of  Zone  1,  q,  as  follows: 
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Where  Allg  is  the  heat  of  decomposition  per  unit  weight  of  binder  plus  the 
enthalpy  change  to  heat  the  materials  in  Zone  1,  fg  is  the  weight  fraction 
of  the  propellant  which  is  binder,  and  is  the  density  of  the  solid 
propellant. 

The  burning  rate  is  expressed  in  terns  of  conduction  and  radiation 
contributions  in  Equation  (1). 


and 


„  4  .P  4^  ,  j  T2  T1  1/3 

r  =  c(T-  -  T,  )  +  d - m  p 

CT-,  +  T,)1/'5 
£  1 


Bbpp  *o 

J  _  1  y,  A  1 

'  2  a0 ' 775 


(6) 

(7) 


Where  Aq  is  the  area  of  propellant  considered  as  a  plane,  the  average 
area  of  propellant  considering  bumps  and  depressions,  e  the  emissivity  of 
the  materials  in  Zone  2,  o  the  Stefan-Boltzman  constant  and  K  the  coefficient 


of  thermal  conductivity  in  Zone  1.  The  constants,  c,  d,  and  T^  were  deter¬ 
mined  from  the  fullscale  motor  burning  rate  data  at  high  pressure,  since 
they  could  not  be  calculated  with  sufficient  accuracy  for  practical  purposes . 
There  is  no  convection  term  since  by  definition  Zone  2  consists  of  materials 
moving  away  from  the  surface.  Further,  the  model  assumes  that  all  time  heat 


transferred  back  through  Zone  2  is  effective  to  pyrolyze  the  binder  and  heat 
the  products  in  Zone  1. 

It  is  important  to  note  that  each  of  the  constants  c  and  d  contain  factors 
related  to  time  surface  area  of  the  propellant.  The  actual  area, including 
bumps  and  depressions,  is  not  necessarily  time  planar  area,  but  is  related  to 
oxidizer  and  aluminum  particle  size  and  concentration,  lime  effect  of  chang¬ 
ing  oxidizer  particle  size,  which  changes  the  surface  area,  is  a  multiplying 
factor  in  the  burning  rate  equation,  and  therefore  time  variations  in  burning 


rate  due  to  particle  size  variations  obsemred  at  high  pressure  can  be  used 
to  calculate  the  variations  at  lom</  pressure. 
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In  using  Equation  1,  it  is  necessary  to  determine  T,  and  T_. 

x  £ 

is  a  constant  equal  to  1200°C  or  about  1500°K,  determined  by  fitting 
the  data.  is  a  function  of  pressure  and  the  expected  pressure  depen¬ 
dence  is  shorn  in  Figure  2.  At  high  pressures,  where  ,  c  and  d  in 

Equation  1  are  determined  from  full-scale  motor  data,  is  the  equilibrium 
propellant  flame  temperature  (calculated  from  the  thermodynamic  considera¬ 
tions)  which  is  lower  than  the  boiling  point  of  aluminum.  Under  these 
conditions,  the  aluminum  boils  and  reacts  rapidly  and  does  not  form  a 
boiling  layer  at  the  top  of  Zone  2.  At  intermediate  pressures,  is  the 
boiling  point  of  aluminum.  At  low  pressures,  where  diffusion  lowers  the 


partial  pressure  of  aluminum,  is  described  by  the  straight  line 
(T  versus  log  P)  tangent  to  the  aluminum  boiling  point  curve.  This  is 
justified  by  a  relationship  between  the  temperature  and  pressure  for  a 
boiling  droplet  of  aluminum  obtained  by  equating  1)  the  rate  of  vaporiza¬ 
tion  by  heat  conduction  to  the  droplet  with  2)  the  rate  of  diffusion  of 
aluminum  into  the  surrounding  bases.  Using  Maxwell's  Diffusion  Law,  the 


relationship  is: 


_  R  iHvT  d£ 

dl  =  •  — p—  • 

F  K  p 


where  R  and  F  are  constants  from  the  diffusion  equation,  AHy  the  heat  of 

vaporization  and  K  the  coefficient  of  thermal  conductivity.  Since  K  is  a 

function  of  T  to  the  first  power  for  a  monotonic  gas  and  AH  for  aluminum 

AH  T  v 

is  virtually  independent  of  the  pressure,  — j-  varies  by  less  than 
5  percent  for  aluminum  between  1500°C  and  3000°C,  and  thus: 

dT  =  or  T  =  A  log  P  +  B  (9) 
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where  A  and  B  are  constants,  The  tangent  point  could  occur  for  some  pro¬ 
pellants  at  pressures  above  the  intersection  of  the  flame  temperature  and 
aluminum  boiling  point  curves ,  However,  the  curves  of  Figure  2  were  drawn 
anticipating  the  results  obtained  below  for  this  specific  propellant. 

Figures  3  and  4  present,  respectively,  the  data  reported  for 
cured,  solid  strands  (approximately  one  inch  diameter)  burning  at  la-: 
pressure,  and  the  calculated  values  for  T2,  using  Equation  10,  where  Tj, 
c  and  d  were  determined  from  the  high  pressure  fullscale  motor  burning 
rate  data  at  the  same  ambient  temperature . 


r  =  2.53  x  10"16  (T24  -  T^)  +  3.20  x  10'4 


T2  ”  T1 


Cr2  +  V 


m p 


1/3 


(10) 


Figure  4  shows  that  at  high  pressure  the  calculated  T?  values  for 
solid  strands  are  400°C  below  the  equilibrium  propellant  flame  temperature, 
and  at  intermediate  pressures  (85  ‘to  45  psi)  are  400°C  below  the  equilibrium 
aluminum  boiling  point.  Further,  the  low  pressure  results  follow  a  straight 
line  on  the  T  versus  log  P  plot. 

These  results  can  be  interpreted  to  mean  that  the  design  of  these 
solid  strands  and  this  particular  solid  strand  burning  apparatus  resulted 
in  inefficiency  in  heat  feedback  relative  to  the  fullscale  motor  which 
lowered  the  heat  flux  to  the  propellant  by  an  amount  equivalent  to  lowering 
T2  by  400°.  Since  the  data  were  reduced,  considering  both  the  radiation 
and  conduction  contributions  tc  the  burning  rate,  this  heat  loss  is  assumed 
to  occur  in  Zone  2. 

Using  T2  calculated  for  solid  strands  plus  400°C,  the  T2  curve  at 
low  pressures  in  Figure  2  was  constructed  as  a  straight  line  tangent  to  the 
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equilibrium  aluminum  boiling  point  curve  at  45  psi.  The  calculated  burn¬ 
ing  rates  shown  in  Figure  5  were  obtained  from  Equation  10  and  Figure  2. 

It  is  seen  that  the  observed  full-scale  data  are  reproduced  exactly. 

Using  the  model  developed  here,  it  is  possible  to  reduce  subscale 
burning  rate  data  where  measurements  are  reported  at  both  high  and  low 
pressures  to  calculate  values  for  T2  which,  when  corrected,  (using  the 
flame  temperature  at  high  pressure  as  the  T0  applicable  to  fullscale  burn¬ 
ing  rates)  yield  predictions  to  the  fullscale  burning  rate  at  low  pressures. 
It  should  be  noted  that  the  high  pressure  burning  rates  and  the 
determined  from  equilibrium  aluminum  boiling  point  information  (assuming 
that  partial  pressure  of  aluminum  is  not  lowered  by  diffusion)  can  be  used 
to  predict  an  upper  theoretical  limit  for  the  low  pressure  burning  rate  of 
the  propellant  in  fullscale  motors  in  the  absence  of  any  subscale  data. 
Further,  the  true  fullscale  burning  rate  is  bracketed  at  low  pressure  by 
the  upper  theoretical  limit  and  any  subscale  data,  thereby  increasing  the 
confidence  in  the  predictions  of  T2  from  subscale  results. 

The  model  developed  here  can  be  used  to  predict  the  lew  pressure 
limit  for  propellant  combustion  for  both  subscale  and  fullscale  tests. 

In  equation  5,  if  T2  =  Tp  there  is  no  exotheimic  reaction  which 
can  feed  heat  back  to  the  surface,  and  therefore  self-sustained  combustion 
is  not  possible  below  the  corresponding  pressure.  The  burning  rate  at 
extinction  is  zero,  however,  at  pressures  just  above  the  extinction  pressure, 
the  burning  rate  can  be  plotted  on  a  log-to-log  scale  as  shown  in  Figures 
3  and  5.  The  extinction  pressures  for  subsonic  flow  are  obtained  by 
extrapolating  the  straight  lines  for  1^  versus  log  pressure  down  to 
T2  =  1500°K,  as  shown  in  Figure  6.  The  lowest  observed  burning  pressures 
are  1.9  psi  (compared  to  1.3  calculated)  for  solid  strands  and  0.8  psi 

(compared  to  0.4  calculated  for  fullscale  motors)  for  a  50  lb  motor. 
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Using  the  pressures  and  burning  rates  determined  here,  and 

calculated  equilibrium  C*  values  for  the  propellant  as  a  function  of 

pressure,  the  throat  areas  necessary  to  extinguish  a  fullscale  motor  when 

supersonic  flow  is  maintained  (C*  extinction)  were  calculated  using  the 

dos 

mass  balance  equation  with  equal  to  zero. 


dp 

v  ~at =  DPV  " 


(11) 


gp 


(12) 


where 

V  =  the  chamber  volume 

pg  -  the  gas  density 
Pp  =  the  solid  propellant  density 
=  the  burning  area 
At  =  the  throat  area 
g  =  the  gravitational  constant 
C*  =  the  characteristic  velocity 

The  values  of  parameters  describing  the  fullscale  motors  were  used 
in  Equation  12. 


TABLE  1 

CALCULATED  THROAT  AREAS 


'ressure 

Burning  Rate 

c* 

Throat  Area 

_(psi) 

(in/sec) 

(ft/sec) 

(sq.  in) 

2.0 

.0165 

4928 

528 

1.0 

.0085 

4888 

540 

0.9 

.0074 

4882 

(522) 

0.8 

.0062 

4876 

(492) 
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The  maximum  value  in  the  calculated  throat  area  occurs  at  about 
1  psi.  At  lower  pressures ,  or  higher  throat  areas,  the  mass  generated  by 
the  propellant  is  less  than  the  mass  leaving  the  chamber  and  therefore 
there  can  be  no  steady-state  pressure  in  the  chamber.  That  is 


dp 


V 


i  I  =  pAbr 


-  PA. 


c£ 


(13) 


for  all  pressures  and  the  pressure  decreases  rapidly  to  the  extinction 
pressure  where  the  propellant  ceases  to  bum.  It  is  noted  that  C*  extinc¬ 
tion  occurs  at  a  higher  pressure  (1  psi)  than  low  pressure  extinction 
(0.4  psi)  in  the  absence  of  supersonic  flow. 


CONCLUSIONS 

From  the  results  of  the  analysis ,  it  is  concluded  that  propellant 
burning  rate  variations  with  pressure  for  aluminized  ammonium  perchlorate 
polybutadiene  propellants  can  be  satisfactorily  represented  by  a  model 
based  on  heat  feedback  rates.  The  description  of  fullscale  motor  propellant 
burning  rates  at  low  pressure  can  be  developed  from 

1.  Fullscale  data  at  high  pressure 

2.  Subscale  data  at  low  pressure 

3.  Physical  chemical  property  data  on  propellant  components 

4.  Thermodynamic  calculations  of  flame  temperature,  C* ,  etc. 

The  approach  leads  to  a  direct  way  of  assessing  the  relative  heat 

feedback  efficiencies  of  subscale  propellant  burning  devices  and  permits 
the  assessment  of  high  pressure  burning  variations  encountered  in 
propellant  production  on  low  pressure  burning  rate.  The  burning  rate 
equation  leads  directly  to  an  accurate  prediction  of  propellant  extinction 
pressures  for  subsonic  and  supersonic  flow  conditions.  Additional  work  is 
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required  to  develop  a  more  fund ament 2 T  derivation  of  the  burning  rate 
relationship  and  to  test  it  on  other  types  of  propellants  with  different 
components,  and  where  exothermic  reactions  occur  in  the  surface,  or  in 
Zone  1. 
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